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Abstract. Sonoluminescence occurs when tiny bubbles filled with noble gas atoms are driven by a 
sound wave. Each cycle of the driving field is accompanied by a collapse phase in which the bubble 
radius decreases rapidly until a short but very strong light flash is emitted. The spectrum of the 
light corresponds to very high temperatures and hints at the presence of a hot plasma core. While 
everyone accepts that the effect is real, the main energy focussing mechanism is highly controversial. 
Here we suggest that the heating of the bubble might be due to a weak but highly inhomogeneous 
electric field as it occurs during rapid bubble deformations [A. Kurcz et al. (submitted)]. It is shown 
that such a field couples the quantised motion of the atoms to their electronic states, thereby resulting 
in very high heating rates. 
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INTRODUCTION 

Sonoluminescence is a phenomenon tliat derives from tlie acoustic cavitation of no- 
ble gas atoms [1]. There are two classes of sonoluminescence: multi-bubble [2, 3] and 
single-bubble sonoluminescence [4, 5]. Single bubble sonoluminescence is character- 
ized by the emission of a strong light flash over a very short period of time from a single 
extremely hot gas bubble. Under appropriate conditions, the acoustic force on a bubble 
can balance against its buoyancy, holding a bubble stable in the liquid by acoustic lev- 
itation. Such a bubble is typically quite small compared to an acoustic wavelength and 
is capable to confine the particles of the trapped van der Waals gas close to their covol- 
ume. For specialized conditions, a single, stable, oscillating gas bubble can be forced 
into such large amplitude pulsations that it produces sonoluminescence during each and 
every acoustic cycle. 

A typical single-bubble sonoluminescence cycle is shown in Fig. 1(a). Most of the 
cycle, the bubble behaves isothermal (c.f. Fig. 1(b)). Point A marks the beginning of the 
collapse phase in which the bubble approaches its minimum radius of about 0.5 /im very 
rapidly with supersonic speed. Here, the bubble becomes thermically isolated from the 
surrounding liquid. In Point B, the temperature within the bubble significantly increases 
with a heating rate of lO^'^ — lO^^K/s and a strong light flash emerges which last for 
about 40 ps. Point C denotes the beginning of the expansion phase in which the bubble 
oscillates around its equilibrium radius until it regains stability. 

The emitted light mainly consist of a continuum of blackbody or Bremsstrahlung 
radiation. Detailed measurements of the light spectra indicate temperatures above lO^K 
[6, 7, 8]. It is even possible to observe light emission in the ultraviolet regime which hints 
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FIGURE 1. A typical single-bubble sonoluminescence cycle, (a): Time dependence of the driving sound 
pressure and the corresponding bubble radius, (b): Time dependence of the temperature. 

at temperatures of about 10^ K in a bubble driven at IMhz [9]. Emission lines from 
transitions between high energy states of noble gas atoms which cannot be populated 
thermally [10, 11] point at the formation of an opaque plasma core [13]. Evidence for a 
plasma core has also been found in multi-bubble sonoluminescence experiments [12]. 

The time dependence of the bubble radius and its nearly adiabatic compression are 
theoretically well understood up to a certain point when it approaches the minimum 
radius [5, 13]. What is the state of the bubble during the last part of the collapse phase 
and the conditions that lead to these enormous heating rates up to very high temperatures 
is still controversial. Here we summarise an idea which suggests that the heating is 
due to the presence of a highly inhomogenous electric field as it occurs during rapid 
bubble deformations [14]. This field couples the motion of the noble atoms to their 
electronic degrees of freedom. When combined with spontaneous emission from the 
atoms, a quantum optical heating process can occur. Similar couplings are responsible 
for the cooling of ions in ion trap experiments [15]. 



Approaching its minimum radius close to point B in Fig. 1(a), the bubble is no longer in 
a thermal equilibrium. Suddenly, an increase in entropy occurs which is based on highly 
irreversible processes. It causes a temperature increase much higher than what can be 
caused by thermodynamic heating processes (c.f. Fig. 1(b)). In the following we address 
two questions: Why does the bubbles need to be filled with noble gas atoms? What is 
the main energy focussing mechanism during the collapse phase of the bubble? 

Close to point B, the mean distance between the noble gas atoms becomes so small 
that interactions between them can be described by a Lennard- Jones potential. Indeed, 
the physical condition of the bubble becomes that of a solid state system. The atoms 
experience an equilibrium between repulsive interatomic forces due to overlapping 
orbitals and attractive forces due to the van der Waals interaction. Thus, any significant 
gain in temperature has to be caused by vibrational motion driven into the quantum 
regime. Furthermore, the presence of light requires the assumption of an open quantum 
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system. 

To model the resulting strong confinement of the atoms, we place each of them into 
an approximately harmonic trapping potential. This allows us to quantise the atomic 
motion during the collapse phase, just before the maximum compression of the bubble. 
Around this point, the motional states of each atom can be described by phonons with 
frequency v. In the next section, we show that the gradient of an electric field inside the 
bubble establishes a coupling between the electronic and the quantised motional states 
of each noble gas atom. The origin of the field can be explained by an inhomogeneous 
charge distribution of ionized species from the dissolved liquid due to rapid bubble 
deformations. 

For simplicity we assume that the atoms are effective two-level systems with ground 
state |0) and excited state 1 1) . The corresponding interaction Hamiltonian contains terms 
that result in the excitation and de-excitation of each atom accompanied by the creation 
and the annihilation of a phonon. Also crucial is the presence of a large spontaneous 
decay rate T of the excited state 1 1 ) which keeps the atoms predominantly in their ground 
state. Although these processes are highly non-resonant, they result in a significant 
change of the mean phonon number per atom and increase the temperature inside the 
bubble by many orders of magnitude, even within a few nanoseconds. 

Suppose an atom is initially in its ground state and possesses exactly m phonons, as 
shown in Fig. 2a. We denote this state by |0,m). Notice that phonons are bosons which 
are described by annihilation operators b with [b,b^] = 1. Consequently, a transition 
into the state | l,m -I- 1) occurs with a rate proportional to y/m+ 1, while the rate for a 
transition into the state 1 1 , m — 1) scales only as ^/m. Since the spontaneous decay rate of 
the atom is relatively large, such a transition is immediately followed by an irreversible 
and predominantly non-radiative transition back into |0). This transfers the atom either 
into its initial state 1 0, m) or into the states 1 0, m — 1 ) and 1 0, m -|- 1 ) , respectively. The net 
effect is an increase of the mean phonon number per atom, i.e. heating, since the phonon 
population in the latter state is higher than the phonon population in |0,m — 1). 

THE TIME EVOLUTION OF THE SYSTEM 

We now consider a single noble gas atom at the position r. This atom is typical for the 
many atoms inside the bubble. Its dipole Hamiltonian equals 

Hint = eD-E(r) (1) 

with e being the charge of a single electron, the (real) atomic dipole moment 

D = DoiC7"+H.c., (2) 

a+ = |1)(0|, <7~ = |0)(1|, and where E is the electric field inside the bubble. For 
simplicity, we assume that all field components point in the direction of a single unit 
vector k. This allows us to write E(r) as 



E(r)=J^E,ei*'- + c.c. 
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FIGURE 2. (a): Level configuration of a single atom-phonon system indicating the immediately rel- 
evant transitions, if the atom is initially in |0,m). £1 and A denote coupling constants and F is the 
spontaneous decay rate of level 1 . (b): The mean phonon number m as a function of time for v — lOMHz 
while a = 10*" Hz, A = 10'2 Hz, r = 10'^ Hz ' , and coo = 10'^ Hz. Good agreement is found between the 
numerical solution of the full rate equations (9) and (10) and Eq. (13) (shaded area). 



with amplitudes E;^ and wave vectors k = kk. Moreover, we consider the atomic motion 
in the k-direction as quantised with b being the corresponding phonon annihilation 
operator. Then k ■ (r — R) = Ax(^b + b^y Here R is the current equilibrium position of 

the noble gas atom with mass M and Ax — ^/hjlMv is the width of its ground state 
wave function in the respective vibrational mode. If the atom is well localized within 
the wavelength of its trapping potential, the Lamb-Dicke approximation allows us to 
assume that exp(ifck- (r — R)) = 1 +ifcAA;(Z7 + Z7^) [15, 16]. Substituting this into Eq. (1), 
we obtain the interaction Hamiltonian 

Hint = m [o- + (7+) + M (Z7 + b^) [a- + (7+) (4) 

with the (real and positive) coupling constants 

n = (2e/^)J^Doi-Re(E,ei*'^«) , 

A = -{le Ax/^) £A:Doi • Im (E^e^*"^ "*) . (5) 

This Hamiltonian is essentially a Jaynes-Cummings Hamiltonian with A being pro- 
portional to the gradient of £l in the direction of the quantised motion of the atom, 
i.e. A = Axk- Vn(R). A strong atom-phonon coupling therefore does not necessarily 
require the presence of a strong electric field. It only requires a highly inhomogeneous 
field inside the bubble. In the following, we neglect interactions between the noble gas 
atoms other than the ones already included in the harmonic trapping potential of each 
particle. Dissipation in form of spontaneous photon emission from the atomic state |1) 
is taken into account by the master equation [16] 
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Here hcoo and fiv are the energy of the atomic state 1 1) and of a single phonon. 

Eq. (6) can now be used to obtain a closed set of rate equations. Its major quantities 
are the phonon number m= (b'^b) and 

^l,2=(CTi,2), X3 = {0+(7--G-a+), Yi = {b + b^), Y2=l{b-b^), 
Y3 = {b^ + b'^^), YA = i{b'-b^^), Zi^2={Oi,2{b + b'^)), Z3,4 = i{Oi,2{b-b'^)) (7) 

with the Pauli operators Ci = (7+ + <7~ and 02 = i(<7~ — (7+). Here we assume 

0)0 > V, r, A and m > 1 (8) 

and approximate the expectation value of operators of the form (Boj,) by (5) ((73). The 
latter applies when the expectation value of B is about the same for an atom in |0) and 
for an atom in 1 1). Eq. (6) then yields 

m = AZ3, X3 = 2(^2X2 + AZ2)-r(X3 + l), Yi = -vY2, 

Y2^2AXi + vYi, Y3 = -2{vY4+AZ3), Y4 = 2{vY3 + AZi) , (9) 

and 

Xi = -cOoX2, X2^-2{a + AYi)X3 + cOoXi, Zi = -cOoZ2, Z3 = 2A-cOoZ4, 
Z2 = -2{aYi + AF3 + 2Am)X3 + cOqZi , Z4 = -2{aY2 + A}4)^3 + (O0Z3 (10) 

up to first order in I/oq. In the beginning of each sonoluminescence cycle the particles 
experience neither a strong trapping potential nor the presence of an inhomogeneous 
electric field inside the bubble. We can therefore assume that the coherences defined in 
Eq. (7) are initially zero and that the atom is in its ground state. Condition (8) allows us 
to simplify the above rate equations via an adiabatic elimination of Eq. (10). Doing so 
we obtain a set of equations where the derivatives of X3, Fi, and Y2 decouple from the 
rest. Solving them for the case of a relatively strong atom-phonon coupling constant A 
with A > n and 4A^ > vcoq yields 

X3(0 = -l, yi(0 = ^^-[cosh(AO-l], Y2{t) = -^-smh{Xt) (11) 

1 /2 

with A = V (4A^/v(Oo — 1) up to first order in I/oq- For times t of the order of 
1/A, Yi and Y2 are of the order of I/cOq. Taking this into account, we find that Zi = 
— 2A(2m + Y3)/cOo and Z3 = — 2AF4/(ao in first order in 1 /(Bq. The variables m, F3, and 
Y4 in Eq. (9) hence evolve according to 

2A2 2(2A2-va>o) ^ SA^ 2(2A2-va>o) 

m= F4, F3 = 74, Y4 = m ^3 . (12) 

(Oo (Oo (Oo (Oo 

For m(0) = mo and 73(0) = 74(0) = 0, this yields 
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m{t) = mo + ' 2 " ^0 sinh^(A?) . (13) 
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As one can see in Fig. 2(b), Eq. (13) describes an approximately exponential heating 
process as long as a relatively large decay F secures that the atom remains in the ground 
state predominantly ^ Taking into account typical experimental parameters, the phonon 
energy in the bubble can easily increase by a factor ten or more, even within a few 
nanoseconds. Using the relation m-hv = ksT, our model can easily predict temperatures 
well above lO^K inside the bubble. 

CONCLUSION 

We attribute the sudden concentration of energy in sonoluminescence experiments to 

the heating of strongly confined noble gas atoms by a highly inhomogeneous electric 
field. The time evolution of each atom is dominated by non-energy conserving processes, 
which result in a permanent increase of its mean phonon number m when combined with 
spontaneous emission. Our model does not contradict current models for the description 
of sonoluminescence experiments, but explains previously controversial aspects of this 
phenomenon. It is based on a quantum optical approach that is routinely used to describe 
the laser cooling of tightly trapped ions [15]. 
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The assumption of a relatively high spontaneous decay rate can be justified by the presence of collective 
effects inside the van der Waals gas formed by the noble gas atoms. 



